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Abstract-The model stability is an important condition, when 
analyzing a financial time series. This assumption may be 
violated because of changes in economical and financial 
conditions. In a parametric model, the model stability 
assumption is reduced to constancy of parameters of model. 
Rolling analysis is useful tool to detect changes in the parameters 
of a statistical model. This paper considers the rolling analysis 
for diffusion processes. Theoretical aspects are proposed and 
some simulation examples are given. A real data set is considered. 
A conclusion section is also given. 
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I. INTRODUCTION 

It is necessary for financial time series to check their 
model stability over time. For example, the volatility of a 
portfolio may increase as well as risk premium rises. This 
concept is also a key assumption for prediction. In a 
parametric model, the model stability condition is reduced to 
model parameters constancy. However, economical 
conditions in real world often change and it isn't reasonable to 
assume that parameters of a model remain unchanged. A 
method to check the constancy of parameters is rolling 
analysis. In this method, over a rolling window of sample, 
parameters estimates are computed. When parameters are 
constant over the sample, these estimates are very close to 
each other. However, if parameters change at some time point, 
then the rolling estimates capture theses instabilities. Some 
excellent references in this field are Colby and Meyers (1988), 
Diebold and Mariano (1995), and Alexander (2001). Zivot 
and Wang (2008) proposed some S+FinMetrics functions for 
rolling analysis in financial time series. 

Although rolling analysis has been received considerable 
attentions in discrete time stochastic process, it seems there is 
no existing literature for continuous time series such as the 
diffusion process. This series is useful tool in many 
applications such as economic, finance and insurance, see 
Kutoyants (2004). Therefore, in the current paper, we 
consider the problem of rolling analysis in diffusion processes. 
The current paper is organized as follows. In section 2, we 
propose theoretical aspects of rolling analysis for a diffusion 
process. Some simulation results and application in a real data 
set are considered in section 3. Conclusions are given in 
section 4. 

II. ROLLING ANALYSIS 

In this section, we prepare the theoretical foundation of 
rolling analysis for diffusion process by presenting some 
simple examples. 

Example 1. Let X t be a Brownian motion with a drift 
defined by 

dX t =/jdt + cjdW t ,te[0,T] 

where X t is standard Brownian motion on [0, T] . This 

model is one of the first candid for modeling the price respectively. 



movement of asset prices. It is clear that we ask if the 
parameters say jU is constant as time passes. Following 
Alexander (2001), a h-period (for some h > 0) simple moving 
average is given by 



A<o>=<l>j>;-.*=^ 



The above quantity is good tool for detecting shift in 
parameter jU . The standard error band around the rolling 

2 

a 

estimate of jU is computed as . A 95% confidence 

h 
interval for jU based on rolling estimate is 



/} t (h)±1.96 






In practice, it is not possible to plot // ( (h)for every t in 

[0, T] . Therefore, we apply this method for some predefined 

time points t j; i = 1, ..., n and we use a discrete approach such 

as Euler method. In application, (7 is unknown and it is 

1 - — 

replaced by its estimate namely (— )/ (X — X) 2 ■ 

n t? '' 

Example 2. Let X t be a geometric Brownian motion 
defined by 

dX t =0X t dt + aX t dW t 

This model is studied by Black and Scholes (1973) for 
option pricing formula. This formula had a great influence on 
financial mathematics. Naturally, changes in #will destroy 
the obtained inferential results if this instability is not detected. 
To check the stability of , it is needed to compute and 
monitor the plot of 
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A discrete approximation for above quantity is 

(-)E ( x *~ Xt ~ ) ■ 

v „ '^— l(i;t-h<t,<t) v v ' 



where X\ h is the number of t t 's belonging to interval 
(t, t + h) and X t is Euler (or any suitable method) 
approximation of X t . The standard error and 95 percent 

(7 a 

confidence band are — = and (h) ±1.96— ;= , 

4h ' Jh 
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Example 3. The option pricing formula of Black-Scholes 
assumes that the volatility parameter is constant. However, a 
stochastic volatility model is fitted much better to financial 
time series. Hull and white (1987) proposed the following 
model 



[dX t =0X t dt + a t X t dW t 
da] = /j<j]dt + ga?dB t 

where W t and B t are two standard Wiener processes such 
that cor(dW t ,dB t ) = p . Let v t (h)= [ -afds . In 

Jt-h 

this case, we should study the behavior of the following 
rolling weighted mean 

h- h x s 
and therefore, we have 

E(9 t (h)) = 9 and var(0 t (/i)) = E{(y t (h)y l ) 

It's easy to see that given V t (h) the 9 Ah) is normal with 
zero mean and conditional variance (v t (h)) .A conditional 
95% confidence interval for 9 , given V t (h) is 



t {h)± 



1.96 

V V rC0 



Example 4. Let X t be an Ito diffusion process indexed by 

parameter 9 satisfying in the following stochastic differential 
equation (SDE) 

dX t =a(X t ;9)dt + b(X t )dW t 

Here, we assume that there is a unique solution for this 
SDE, see Kutoyants (2004). Let h be the maximum 
likelihood of based on observations X s , S ^[t — h,t] . 
Following Kutoyants (2004), as h — > co ; the estimator 
0. h is consistent estimator for and it is also unbiased. 

Therefore, we can use th as a shift detector. Also, this 

estimator is asymptotically normal for each fixed t, therefore, 
we can construct an asymptotically confidence bands. If 
was a vector of unknown parameters again, the above method 
works. For example, in Black-Scholes model (example 2), 
using the rolling maximum likelihood of (0, <x) we can 
detect the changes in two parameters, simultaneously. 

Example 5. In example 4, suppose that the diffusion 
coefficient <7 may also relate to unknown parameter In 

this case, it is enough to apply the Ait-Sahalia (2002) 
transform 



Y t =y(X t ;0) = \ 



x, du 



b(u;0) 



To derive a SDE with unit diffusion coefficients in the 
form of 



dY t =a Y (Y;0)dt + dW t 

and again we apply the approach of example 4. Here, 

b( r -'(y,ey,f>) 2dx 

For example, consider the SDE 



dX t =0X t dt + z X t dW t 



rir,0) = \ 



y du 1 w , a(x;0) _ 1 



-My), 



and 



b(x;0) 



a(x;0) = . Therefore, a ( y; Q) = 'LJL and 

dx 



dY, 



2 
applicable here. 



dt + dW ■ Again, approach of example 1 is 



III. SIMULATIONS 



In this section, we propose some simulation results and a 
real data set are studied. 

A. Simulations. 

In this sub-section, we present some simulation results and 
a real data set analysis to evaluate the performance of rolling 
method appeared in section 2. We refer to simulation schemes 
in Zivot and Wang (2008). We consider two cases, (a) the 
Black-Scholes model under the stochastic volatility 
assumption and (b) monitoring the stability of drift and 
diffusion coefficients of a SDE, simultaneously. 

Case (a). Consider example 2, 

dX t =0X t dt + aX t dW t ,O<t<l. 

First, we let X n =\0=O,CT=\t= ,Z =],..., 400 and 

400 

x t . - x t . 

n h = 24 We also compute returns V i = — '■ '—. The 

return series r ; is plotted in Figure 1 (page 7). The standard 
error of t (h) and <J, (h) are approximated by 

t (h) a t (h) 

'. — and ' , 

yj n h 4 2n h 

respectively. Figure 2 (page 7) shows the rolling estimates of 

0, (h) with 95% confidence intervals. Note that since we are 

monitoring and <J simultaneously, to avoid excess of 
error, according to the Bonferroni law, we obtain the marginal 
intervals for each parameter in significant level 975%. The 
rolling drifts and diffusions clearly vary over the sample. The 
rolling drifts start out around 0.1%, fall close to 0% very soon, 
rise again a little, but we can assume that it is zero for the 

entire sample. The plot of C t (h) is also derived (which isn't 

given here). The rolling diffusion values start out around 1 
and remain unchanged for all sample points. These facts 
motivate us to assume that the parameters are constant over 
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the sample. Now, suppose that a change has occurred in 
and it equals to until t 200 and shifts to 0.7 after f 200 Then 

the plot of t (h) (Figure 3, page 7) obviously shows that 
there is a shift in parameter 6 




■:: 



Fig.l Return of B-S series: no change 
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Figure2. Rolling means with confidence bands 
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Fig.3 Rolling means with a change point 



Case (b). In example 3, we derived the conditional 
confidence interval using rolling drift estimates. However, in 
real situations, observations sampled from the stochastic 
volatility are not fixed. To avoid this difficulty, we should 
derive the quantile of marginal distribution of 

t (h),te[O,l\Le.,F(z) 

F(z) = E{®(j^h)(z-0))}, 

where 0(.) is the cdf of standard normal distribution. For 
example, for 6 = 0.3 , jU = 0.4 , £ = 1, for t = 0.5 and 
h = 0.05 , Monte Carlo estimate of F{z) with M = 1000 
repetitions is 0.45. Related standard deviation of Monte Carlo 
estimate is 0.0157. The 0.95-th quantile of t (h) is the 
solution of 



E{0(V^)(q o . 95 -#))} = 0.05 
A Newton-Raphson method is applied here as 



k _ ft— 1 

"o.95 — Wo.95 



E{fiXh)0(JT(hj(q k o ; 9 l-d))} 



where 0(.) is the density of standard normal distribution. 
Using a proper numerical software, we found that Q = 3.5 . 

*0.95 

B. Real Data Set 

The diffusion process has many applications in modeling 
the economical and financial time series. In this way, Hansen 
(1992) considered the annual U.S. output growth rates over 
time period 1889- 1987. He fitted a discrete time AR(1) 
model for this series. He examined many goodness of fit tests 
and concluded that his model worked very well for the data 
set in hand. He also derived many prediction and concluded 
that again this model works very well. He showed that this 
model has remained stable. Therefore, there is no change 
point in this model. However, since every discrete time series 
is an approximation of a continuous time series, following 
Hansen (1992), we also fitted an continuous time AR(1) (as 
alternative and approximation method) to this series in the 
form of 

dX t =aX t dt + adW t . 

We found that a = -0.68 and a = 0.045 . We applied 
the above mentioned method for this data set and we 
understood that the maximum of absolute of difference of 
rolling drifts and CC is less than 0.1. This shows that the 
parameter CC remains unchanged during period of study. 
Therefore, since the continuous time series has remained 
unchanged, the discrete time series has also no change point. 
We also proposed many goodness of fit tests, where their 
results were not given in this paper. 

IV. CONCLUSION 

Diffusion processes play important role in statistical 
inference and the model stability is a necessary assumption 
for them. Rolling estimates monitor the constancy of 
parameters of a model over time. In this paper, we considered 
the rolling analysis in a diffusion process. Theoretical aspects 
are given and by some simulations, we show that these 
estimates detect changes very well. We proposed some 
examples and a real data set. Our real data set is Hansen's data 
set about annual U.S. output growth rates over time period 
1889- 1987. We also understood that the series remains 
unchanged during the mentioned period of time. 
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